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The switching of magnetic layers is studied under the action of a spin current in a ferromagnetic 
metal/non-magnetic metal/ferromagnetic metal spin valve. We find that the main contribution 
to the switching comes from the non-equilibrium exchange interaction between the ferromagnetic 
layers. This interaction defines the magnetic configuration of the layers with minimum energy and 
establishes the threshold for a critical switching current. Depending on the direction of the critical 
current, the interaction changes sign and a given magnetic configuration becomes unstable. To model 
the time dependence of the switching process, we derive a set of coupled Landau-Lifshitz equations 
for the ferromagnetic layers. Higher order terms in the non-equilibrium exchange coupling allow the 
system to evolve to its steady-state configuration. 

PACS numbers: 75.70.-i 75.60.Ej, 72.15.Gd 75.10.-b 



I. INTRODUCTION 

The possibility of using the exchange field of a spin- 
polarized current to aid in switching of a magnetic layer 
is not only an intruiging prospect for future applications 
in the field of magneto-electronics but also a challenge 
with interesting physics. There is still a lot of ambiguity 
in the interpretation of the complex data of the recently 
observed switching of domains in spin valves j]J and tun- 
nel junctions || , of magnetic clusters || an d layers j| , 
nonetheless, in Refs. Jl],|3|,f|] it has been argued that the 
switching occurs through relaxation of conduction elec- 
tron spin-polarization to the local moments of the ferro- 
magnetic layers as proposed by Slonczewski Mffl. 

In this work we introduce a different model for the 
switching, where the spin-polarization does not relax. 
The effect is that the spin current carries an exchange 
field that acts on the local moments of a layer and forces 
its magnetization to take the orientation of the spin- 
polarization of the conduction electrons. 

In the conventional view one considers electrons flow- 
ing in the direction of the net current from, say, a "fixed" 
layer that has a large magnetic moment to a "free" layer 
that has a small magnetic moment which is easy to re- 
orient (see Fig. [j]). The effect on the free layer is indeed 
relatively small as the associated energy, which is of the 
Zeeman-type, is proportional to the product of the small 
moment of the free layer and the exchange field of the 
conduction electrons polarized by the fixed layer (refer 
to upper part of Fig. |). However, this picture neglects 
the much larger, albeit counter intuitive, effect of spin- 



dependent reflection of electrons by the free layer, so that 
spin-polarized electrons move in the direction opposite to 
the net current and interact with the large magnetic mo- 
ment of the fixed layer. The energy of this contribution 
is thus also much larger and of opposite sign (see lower 
part of Fig. |). 

Taken together, both contributions establish a strong 
non-equilibrium exchange interaction (NEXI) between 
layers 0. The sign of this interaction is determined by 
the direction of the current; it does not oscillate and its 
range is controlled by the spin diffusion length of the con- 
duction electrons |8||| . Therefore, the NEXI is a volume 
effect and exerts a torque (force) throughout the volume 
of the magnetic layers (leading to precession); it is the 
dominant mechanism that drives the system towards its 
switching threshold. Here, we describe a simple model 
that shows how the NEXI defines the magnetic config- 
uration of the layers with minimum energy, while the 
relaxation of the conduction electron spin polarization 
provides a way for the system to evolve to its minimum 
energy configuration. 

The rest of the paper is structured as follows. In Sec. II 
we introduce our model of a spin valve and give a brief 
derivation of the NEXI. To illustrate the switching mech- 
anism of the spin valve, in Sec. Ill we derive a set of 
coupled Landau-Lifshitz equations for uniform magnetic 
dynamics, whose structure is similar to those of antifcr- 
romagnets or ferrimagnets; however with the important 
difference that they are current driven and the magnetic 
sublattices are replaced by the spatially separate mag- 
netic layers. The following stability analysis in Sec. IV 
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shows that depending on the direction of the current ei- 
ther the parallel or the antiparallel configuration becomes 
unstable beyond a critical current, so that, in principle, 
even an arbitrarily small amount of relaxation allows a 
switching of the magnetically "softer" , i.e. free, layer. De- 
spite the simplicity of our model, quantitative estimates 
are in reasonable agreement with experiments. In Sec. V 
we introduce Gilbert damping to describe possible re- 
laxation processes that allow the system to evolve over 
time to its minimum energy configuration, i.e. to switch. 
Finally, we compare our results to the models by Slon- 
czewski |3]||, and Berger in Sec. VI and conclude in 
Sec. VII. 
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FIG. 1. The geometry of a standard trilayer spin- valve; two 
planar ferromagnetic metal layers of thickness Il(r) and total 
magnetic moments M i ( J j) at an angle <d relative to each other 
are separated by a non-magnetic metal spacer NM of thickness 
d. 



II. NON-EQUILIBRIUM EXCHANGE 
INTERACTION (NEXI) 

The geometry of our model system is a typical spin- 
valve structure shown in Fig. [l]. It consists of two planar 
ferromagnetic metal layers of thickness Il(r) whose total 
magnetic moments M^m are at an angle relative to 
each other. The magnetic layers are separated by a non- 
magnetic metallic spacer NM of thickness d. A complete 
treatment of the magnetic dynamics of such a spin valve 
with a current in the perpendicular direction requires a 
simultaneous solution of the equation of motion for in- 
dividual spins of the magnetic layers (see for example 



Ref. ) , and equation of motion of the spin-polarized 
charge carriers B| . 

To focus on the essential effects, the calculations are 
made on the simplest assumptions; we model the poly- 
crystalline thin films as uniformly magnetized layers with 
uniaxial anisotropy in the plane of the layers. For clarity, 
we do not focus on the details that establish the equi- 
librium (zero current) coupling between magnetic lay- 
ers; while the experiments Q suggest that the equilib- 
rium Ruderman-Kittel-Kasuya-Yosida (RKKY) coupling 
is negligible, the omission of the dipolar coupling between 
the layers (fringing fields) is probably a gross simplifica- 
tion of the physical picture. We treat the steady state 
non-equilibrium (constant current) situation in terms of 
the NEXI which can be written as a sum of quantum- 
interference and current driven terms; see for example 
Eqs. (2) and (3) in Ref. §. Although the first contribu- 
tion to the coupling is finite at equilibrium, a.k.a. the 
RKKY interaction, it is a "surface" effect because it os- 
cillates and scales as 1/d 2 [OJ, therefore we will neglect 
it here as pointed out above. We posit that it is only 
the current driven term, which is a volume effect, that is 
responsible for the switching of the layers. Its decay is 
controlled by the spin diffusion length which is consider- 
ably longer than the typical thickness of spacer layers in 
metallic multilayers. Although our calculations are taken 
in the ballistic limit, they can be generalized to account 
for diffuse transport as will be pointed out in the text. 

When a constant current is driven in the direction per- 
pendicular to the plane of the magnetic layers, it becomes 
spin polarized. We define a spin current as: 

*U%r))=3%r)*HR), (i) 

whose polarization is along the unit vector of magneti- 
zation n-L(R) — m L(R)/ m L(BA^°£ the local moments, i.e. 
along the z^^-axis in Fig. [y, which generates the po- 
larization, and whose magnitude is 

3l(r) = — {JL(r) -r L (R)) = — -?c Vl(r) ■ (2) 

We defined the electric charge and Bohr magneton as 
e = — |e| and \ir. = \e\h/(2m*c), respectively; m* being 
the effective mass of the conduction electrons. Outside a 
magnetic layer the spin current decays within a distance 
of the spin-diffusion length A s f , so that t]l(r) is propor- 
tional to exp(— a;/ A s f) where x is the distance away from 
the layer. The factor t]l(r.) describes the spin-dependent 
reflection (for diffuse transport this leads to the effect of 
spin accumulation and can be included in rjL(R)) and can 
vary between -1 and 1. 

The NEXI comes then from the coupling of the spin 
current generated by the left layer S(j^) interacting 
with the local moments in the right layer and vice versa. 
When we take as the direction of positive current from 
left to right, the coupling in linear response is: 

El* xi = E R (jM)+E L (-j%). (3) 
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The local non-equilibrium coupling El(r) is proportional 
to the scalar product in spin space of the spin current and 
the local moments, i.e. Er oc ^s(jff) ■ and similarly 
for El, where we averaged the spin-operators over the 
non-equilibrium statistical ensemble of the entire system, 
and have neglected spin-fluctuations [pLlJl • From Eq. (|l|) 
follows, 3f(-i*[ H) ) = so that E L {-j%) = 

—El(Jr), and Eq. (||) takes the form: 
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(4) 
are the 



local non-equilibrium exchange fields in the mean field 
approximation. Here Jl(r) is the coupling constant be- 
tween conduction electrons and local moments; Hl(r) = 
9l(R)I j 'Bi and guR) is the Lande factor in the respec- 
tive layer. The preceding discussion can be visualized 
by Fig. 0. Another way of writing Eq. (|J) is in terms 
of an effective Heisenberg coupling that depends on the 
current: 
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M R . (5) 

In Eqs. (|i|) and (|J), respectively, we also introduced 
the simplifying assumption that the amount of spin po- 
larization 7?l(_r) does not depend on the angle O between 
the magnetic layers (see Fig. [I]) and that there is no relax- 
ation other than spin-diffusion. In fact the spin-polarized 
current of each layer depends on boundary condition to 
adjacent layers, e.g., as outlined in Ref. Q, so that 
Vl(r) h as to be calculated self-consistently from Eqs. (||) 
and (Q), by taking into account the back effect of the cou- 
pling between the layers on the spin-polarization itself. 
A simple illustration of an approach which is self con- 
sistent in the current but not in energy can be given in 
the limit of diffusive transport by calculating the differ- 
ence in energies for parallel and antiparallcl configuration 
associated with the spin accumulation in the spin-valve. 
Using the picture developed in Ref. JI3], one finds that 
for different layer thickness this energy contribution is 
non- vanishing. Details will be given elsewhere for a fully 
self-consistent calculation. 

The form of Eq. (Q) leads to some immediate conse- 
quences. If the system is symmetric in its magnetic prop- 
erties, then, for the assumed linear response regime of 
the current, which is reasonably well satisfied in metallic 
multilayers even for high current densities, the two cou- 
plings El and Er are equal and opposite, so that there 
is no non-equilibrium coupling between the layers. If, 
however, the magnetic properties or the thickness of the 
magnetic layers are dissimilar (in general, this would in- 
clude non-uniform magnetization), the spatial symmetry 
of the system is broken and a current driven coupling 
between the layers exists. In particular, on reversing the 
current the interlayer coupling changes sign. 
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FIG. 2. The figure shows the two contributions of the 
spin-current to the non-equilibrium exchange interaction 
(NEXI) given by the hatched regions: the forward contri- 
bution jjf from the left ferromagnet and the reflected contri- 
bution from the right ferromagnet. The NEXI between 
the two ferromagnetic layers, as shown here for parallel con- 
figuration, is given by the difference of the double hatched 
regions which are proportional to jff Mr and j'Jf Ml, respec- 
tively. Only in certain limiting cases one of the contributions 
can be neglected. 



III. EQUATIONS OF MOTION 

To obtain the equation of motion for the magnetization 
of the spin valve, we follow the standard procedure and 
define an effective field (ll|, 



ueff _ h (i) 
L(R) - L(R) 



^L(R)^L(R) ( m L(R) ■ n L(R) 



(6) 



which is derived from the total energy of the system in- 
cluding the effects of uniaxial anisotropy, dipole-dipole 
interactions, and -E^exi °f Fq. (Q) . The internal field (first 
term in Eqs. (^J)) includes a contribution from the exter- 
nal field h.( c \ the spin current, and the magnetic dipole- 
dipole interaction which in a uniformly magnetized ellip- 
soid takes the form 
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L 
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h (e) + h noxi _ 4^ . mL ( 7a ) 
h (e) + h ncxi _ in ^ R . mR ( ?b ) 

where Nl(r,) is the tensor of demagnetization factors, and 
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are the effective non-equilibrium coupling fields between 
the layers induced by the spin-current, which vanish for 
identical layers. In general, when magnetic layers cannot 
be described by single domain particles, this symmetry 
condition does not hold due to the non-uniformity of the 
layers. Further, in such a case the coupling fields (||) will 
contain additional contributions from the fluctuation of 
the magnetization. The second term in Eqs. (||) is due to 
the uniaxial anisotropy, where Pl(r) is the corresponding 
constant and its in-plane unit vector. We would 

like to point out that the assumption of uniform magneti- 
zation means, in addition, that the effect of the induced 
magnetic field from the current, which leads to vortex 
formation, is small compared to the coupling field h noxl 
so that we can use the magnetostatic approximation in 
which h( e ' is the field produced by currents in coils alone. 
The equations of motion are then given by 
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dt 
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(9) 



where Jl(R) = 9l(r)H>b /% > are the gyro-magnetic 
ratios. The damping terms R-l(R) wn l De discussed in 
detail later. 

Parenthetically, the solution of Eqs. (^) is closely re- 
lated to solving the equations of motion for an antifcr- 
romagnet or ferrimagnet, where now an effective field of 
the exchange coupling does not act between sublattices 
but different layers. Therefore, due to the NEXI there 
exist long wavelength "acoustical" and "optical" modes 
wherein the magnetizations precess either in or out of 
phase. In an experimental set-up similar to the one de- 
scribed in Refs. |l6|-|T8f|, this would offer a direct pos- 
sibility to measure the coupling between the magnetic 
layers as a function of applied current. The details will 
be presented elsewhere. 



IV. CRITICAL CURRENT 

So far we have derived the equations of motion for two 
coupled magnetic particles. The solution of the equa- 
tions requires to divide the problem into a stationary one, 
which we shall consider first, and time-dependent one. 
In other words, before studying any form of dynamic be- 
haviour of the magnetization, one must first determine 
the steady-state, i.e., the orientations of the magneti- 
zations in the absence of time dependent driving fields. 
This orientation depends on the current. In combination 
with Eqs. @, we obtain the equilibrium (steady state or 
constant current) configuration of the moments from the 
conditions: 
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^L(R)\ 
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0, 



(10) 



which are independent of the time-dependent damping 
term and where also the effective field h^f R * | q does not 
depend on time. These conditions have to be satisfied 



simultaneously for both magnetic layers. For many ap- 
plications, there exists an interesting limiting case where 
the thickness and anisotropy of one of the layers is very 
large, termed the fixed layer, so that one can assume, for 
example, = n^. In other words, the back effect of the 
right layer, termed the free layer, will be negligible on 
the orientation of the left layer but not on their mutual 
orientation. To be more explicit, we study the steady 
state of this example in more detail. 

Instead of applying the conditions in Eqs. (|l0|), a more 
straightforward method is to minimize the total energy of 
the system. By using the notation we introduced in Fig. |l| 
and assuming that the geometry of the film constrains 
the magnetization to be in the plane of the layer, the 
total energy of the system, in the absence of an external 
magnetic field, depends only on the angle between the 
magnetic moments of the left and right layer: 



0. 



E{Q) = - "-fm L M L - Zfm R M R cos^(e - Q%) 
-h R exi M R cos6, 



(11) 



where @|j is the angle between the easy axis of the right 
layer with zl, and 
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(12) 



is the effective field value of the NEXI on the right layer of 
Eq. (8b) with (1l = V>R = A 1 - From Eq. ( |l2"|) it follows that 
for asymmetric magnetic layers a coupling exists which 
will be dominated by the spin current generated by the 
right layer if the total magnetic moment of the left layer 
exceeds that of the right one, i.e. Ml > Mr. In recent 
experiments one of the magnetic layers is chosen to be 
much thicker jjj so that Ml 3> Mr, and in the following 
discussion we shall assume this is the case. From the 
discussion after Eq.(l) it should be noted that h s £^, RL ^ 

is proportional to X s {/Ir(l) f° r ^L(R) > -W inasmuch as 
jljR) 15 limited by the spin diffusion length A s f. 

It is insightful to express the proportionality between 
the field h R cxl in Eq. (pj) and the current density which 
generated it by the following relation: 



^ ncxi 
Anexi _ ,l R 

A R = —— 

Je 



1 / JlVr M l _ JrVl 
Lie\v[ Mr v% 
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(13a) 
(13b) 



Eq. (13b) holds if we assume that the polarizations are 



the Pauli factors rj L ( R ) = >/L(fl)™L(i?.)/( 4e F ^ M)j and 
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as: 



the Fermi energies. One can, then, rewrite Eq. (HJ 



incxi _ a ncxi ■ 
n R — IV R Je- 



(14) 



In other words, h R exl plays the role of an exchange bi- 
asing field on the right layer generated by the current 
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through the system. To get a better understanding, we 
give a rough estimate to the magnitude of A™ xl . If we 
assume a magnetization for Co of run = 1-42 x 10 6 A/m, 
a Fermi velocity of about 1.5xl0 6 m/s, the layer thick- 
ness rat io (I I — Ir)/Ir=3, and use the free electron mass 
in Eq. ( 13b ), the only undetermined parameter is ??l(_r)- 
Taking the hypothetical case that the amount of spin- 
polarization is the same throughout the system as in the 
ferromagnetic layer, we estimate for Co with r/ = 38% 
A^P X1 to be approximately 5xl0 -5 m, which amounts to 
l/A n R cxi 0.02 mA/kOe. The latter is about an order 
of magnitude lower than the value measured in the ex- 
periments on Co/Cu/Co pillars g. Therefore, our esti- 
mate for A^ cxl derived for perfect spin polarization can 
be seen as an upper threshold for the proportionality be- 
tween biasing field and current. Realistic estimates are 
obtained by adjusting the value for r/ L ^ taking spin- 
diffusion, spin-dependent reflection at the interfaces, and 
the resistance of the layers into account f7],[l4[ . Then for 
fj = 10% [ fL9| , we obtain direct agreement with experi- 
ments l/A£P xi « 0.29 mA/kOe g. 

Having introduced the effective field h * R X1 of the NEXI, 



we realize from the form of Eq. ( Jl 1| ) that finding the 
switching threshold of the right layer (free layer) can be 
treated as if the layer were a single magnetic domain in an 
external magnetic field. The equilibrium direction of the 
magnetization of the right layer (free layer) is then ob- 
tained by the extremum energy condition (dE/d&) = 0. 
In addition, for the equilibrium to be stable (unstable) 
the following relation must hold: (d 2 E/d6 2 ) > (<) 0. 
At (d 2 E/d® 2 ) — the model predicts a transition from 
a gradual rotation of to a sudden switching towards 
the direction of n^, i.e. an irreversible magnetization ro- 
tation, which determines the critical coupling field /i" exl . 
Since we have two conditions and the two unknown vari- 
ables and /i" cxl , we can eliminate and derive an ex- 
pression for the critical coupling field /t" cxl obtained from 
the following relation (refer to Ref. pOl for the details of 
the derivation): 



sin 26^= — 

' V 1 



1 /4-y 2 



3/2 



(15) 



where y = — 2h^ exl /(PRm R ). From Eq. ( |l5| ) we find 
that the critical field is maximum at O^ = and tt/2, 
where, /i" OXI = —f3nm R , and minimum at Q R = 7r/4 
when /i" oxl = —j3nm R /2. For a parallel orientation of 
the uniaxial anisotropies between the layers <d R = this 
translates together with Eqs. (^) and ( |l2| ) to the following 
condition for a spin current induced switching: 
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A R m R 
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m* t]LT]R l L vi - Irv r 
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(16a) 
(16b) 
(16c) 



where j% is the critical current density, and K R = 
Pr m j?/2 the standard expr essio n for the uniaxial 
anisotropy constant pi]. Eq. ( |l6b| ) holds if we assume 
again that the polarizations are the Pauli factors, and 



Eq. (16c) is applicable when Ml ~> Mr. From Eq. (16c) 
it becomes clear that the switching threshold is deter- 
mined by the coupling of the spin current generated in 
the right layer (free layer) to the magnetization in the left 
layer (fixed layer); it would be inappropriate to replace 
the left layer in the problem by a simple spin-polarized 
current source in th e swi tchin g experiments jlj]. On the 
contrary, Eqs. (16a) and ( |l6b| ) show that for almost iden- 
tical layers the critical current becomes very large as the 
denominator tends to zero, so that only transient pro- 
cesses are relevant and, indeed, switching has not been 
observed |22|. Since r\h^R is a measure of the coupling 



strength, the critical current j% is thus inversely propor- 
tional to the coupling of the free layer to the fixed layer, 
so that a strong coupling reduces 011 the other hand, a 
high value of anisotropy K R will increase More gener- 
ally, is proportional not only to the uniaxial anisotropy 
but to whatever constrains the local moments; for exam- 
ple jg will also be strongly influenced by the dipolar cou- 
pling between layers and quantum interference (RKKY) 
contributions to the interlayer exchange coupling. Al- 
though K R certainly cannot capture the details of the 
switching, it yields a simple analytical solution that pro- 
vides the intuitive picture that current driven switching 
can be treated as a single magnetic particle in the current 
driven exchange field h R cxl . This description is particu- 
larly appealing as it allows one to treat generalizations of 
our model similar to those known from magnetic record- 
ing f§. 

We chose = to be the parallel alignment of the 
magnetic moments for zero current and the direction of 
current from left to right which led to a negative sign 
in Eq. (|TT1). Thus, the parallel con- 
h R exl is positive. Having 
to the current (refer to Eqs. (||) and (fL2|)), 
the system becomes then unstable for a sufficient nega- 
tive current ( |l6|) which forces the spin valve to switch to 
an antiparallel configuration. Similarly, starting from an 
antiparallel orientation, a positive current switches the 
spin valve to parallel. We can also give a rough estimate 
to the magnitude of required critical current densities. 
Using again the same data as before of Ref. M and an 
uniaxial anisotropy for Co of K R = 1.5 x 10 4 J/m 3 , we 
estimate the lower limit of t he cr itical current for per- 
fect spin polarization in Eq. (|l6b|) to be approximately 



in front of h R cxl 

figuration is only preferable if h R exl is positive 
related h n R cxi 



0.4xl0 7 A/cm 2 ; this is more than an order of magnitude 
lower than what is needed for the actual switching ob- 
served in the experiments Q , consistent with our limiting 
estimate on 1/A^ CX1 . Using, however, the more realistic 
value ?y = 10% we obtain a critical current of 0.6xl0 8 
A/cm 2 in good agreement with the experiments B. 
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V. RELAXATION PROCESSES 

The conditions for a critical current were obtained from 
the solution of the stationary problem. To describe the 
dynamics of the switching process, one has to solve the 
time-dependent problem. In addition, we have to intro- 
duce relaxation, because without it the right layer would 
only precess but never switch, no matter how strong the 
current would be. However, when the steady state config- 
uration is unstable, even the slightest relaxation is suffi- 
cient for the magnetization to switch. A wealth of possi- 
ble relaxation mechanisms arises from the non-uniform 
motion of the magnetic moments, which we have ex- 
cluded here by making the assumption that the layers 
can be treated as two coupled uniformly magnetized par- 
ticles; also we neglected the back effect of the NEXI on 
the spin polarization ri L ^ R ) which reduces the coupling. 
An important mechanism for relaxation is the transfer 
of angular momentum from the conduction electrons to 
the local moments as discussed by Slonczewski and 
Berger 0. 

Although it is possible to derive this contribution by 
systematically taking the perturbation expansion of the 
s-d Hamiltonian to third order, here we will simply posit 
it by introducing a phcnomcnological Gilbert damping in 
Eqs.®: 
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dM 
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L(R) 



at 



(17) 



where oll(r) is a dimensionless damping parameter that 
depends on the current and on the layer thickness. We 
concentrate only on the part of ffM.ug\/&t due to the 
NEXI; all other contributions follow in analogy from the 
expression of the effective field (roh. If we substitute in 



Eqs. 



) 7i? with 7^ = 7_r/(1 + Or) we can transform 



Eq. ([17 ) into one of the following forms: 
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MrX (Mr X M L ) 

= -ujr (m R cos 8 - x^h'iT) V R , 



(18a) 
(18b) 



where lor — aRj R h R 0X1 is the damping frequency due 
to the NEXI, Xr — m R/h R cxl the respective susceptibil- 
ity, and Vr the volume of the layer. The first way of 
writing the relaxation resembles the effect of the conduc- 
tion electron spin relaxation of Refs. (|J(|. The result 
therein, that the surface torques of both layers impel 
and m.R to rotate in the same direction in the plane of 
the layer, applies also to our case as can be seen from 
Eqs. @ and (0). However, our conclusions are differ- 
ent. With reference to Eq. (^|), we note that the NEXI, 
being the leading order contribution to the switching, 
determines the precession so that the effect of is 
to reduce the coupling and thus minimize the energy in 
Eq. ([Tl]). This is demonstrated by the second way of 
writing Eq. (|l^) which leads to the interpretation that 
niR relaxes towards the non-equilibrium exchange field 



h.R CX1 whose orientation is determined by that of the con- 
duction electron spin polarization which is controlled by 
rri£ and the direction of the current. In other words, if 
the system becomes unstable, R^ 0X1 dominates all other 
contribution in the Gilbert damping, so that R™ xl drives 
the system to its new minimum energy configuration, i.e. 
towards the direction of h??* 1 . This reasoning is analo- 
gous to what happens to a single magnetic particle in an 
externally applied switching field. 

Finally, the way in which we introduced the phe- 
nomenological damping, relates the time for the total 
magnetization to reorient to the relaxation frequency lor; 
the switching motion becomes more viscous as ojr be- 
comes large. This coincides with the common notion that 
the switching time is expected to be fastest for a mod- 
erate value of ujr, and might allow one to use Eqs. (|^) 
despite their simplicity to model current driven switching 
dynamics in more realistic systems. 



VI. DISCUSSION 

We have shown that in order to describe current driven 
switching of magnetic layers in sub-micron sized spin- 
valve structures, it is in general necessary to include the 
non-equilibrium coupling between layers. Only in cases 
where the system is in linear response and the mag- 
netic layers are identical, or the spin-diffusion length very 
short, one can neglect its contribution. The assumptions 
taken in our model are that of two single domain parti- 
cles with an uniaxial anisotropy and a non-equilibrium 
coupling mediated by the spin-current and a local s-d 
exchange interaction. These are certainly gross simpli- 
fications of the physical situation, however, they give 
quantitative agreement with experiments for a reason- 
able choice of the magnetic parameters. 

In making comparison with Slonczewski's model |&^|, 
it should be noted that we make a different set of as- 
sumptions. In his case the NEXI was neglected and the 
system was reduced to consisting of a magnetic layer, 
that serves as a constant source of spin-polarized elec- 
trons, and a single domain particle which acts as a per- 
fect spin filter, i.e. that all majority spins are transmitted 
whereas all minority spins are reflected. The assumption 
of spin-filtering distinguishes Slonczewski's theory from 
that of Berger |l(J , who invokes inelastic spin-flip scatter- 
ing. However, only inelastic spin-flip scattering allows for 
multiple spin- wave excitations that can reorient the axis 
of quantization in a magnetic layer and serves as a dif- 
ferent type of switching mechanism in the limit of short 
spin-diffusion length, as will be shown further down in 
the text. 

Spin-filtering, which we termed spin-dependent reflec- 
tion, gives rise to interlayer coupling. It is non-dissipative 
to the leading order in the coupling which is described by 
the NEXI. Only the next order term in the interlayer cou- 
pling is proportional to Eqs. (p"8}). By neglecting the term 
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in le adin g order, Slonczewski finds only terms similar to 
Eq. ( 18a ) as the sole origin of multiple spin- wave excita- 
tions as predicted by Berger |p3| . Our model produces 
the opposite effect: the precession of is damped and 
relaxes towards the non-equilibrium exchange field h™ XI . 
It is also not surprising that after fitting the phcnomeno- 
logical damping parameter ac |0 ano - adjusting rj to 
10% JkJ , which coincides with our estimate, a good quan- 
titative agreement with experiments is reached; the com- 
petition between the different damping terms Rr just 
reflects the competition between the different effective 
field contributions h^ F (||) in the Landau-Lifshitz equa- 
tion (|). 

Nevertheless, Slonczewski's theory is insightful in that 
it can describe spin-wave instabilities and shows how 
switching can occur due to inelastic spin-flip scatter- 
ing at a multilayer interface after adaptation accord- 
ing to Berger p0| . To understand how such a switch- 
ing mechanism is in principle possible, one has to con- 
sider first the creation of spin-waves in magnetic mul- 
tilayers which were initially observed in point contact 
measurements of Co/Cu multilayers above a critical cur- 
rent As pointed out before, at the root of this 
phenomena lies a relaxation mechanism between conduc- 
tion electrons and local moments due to spin-flip scat- 
tering. When a current flows from a non-magnetic metal 
layer to a ferromagnetic one, its distribution over spin- 
up and spin-down currents has to change. Given that 
most scattering events will conserve the current in each 
spin-subband, a difference in "chemical potentials" ap- 
pears on the scale of the spin-diffusion length away from 
the interfaces |25) which leads to the effect of spin accu- 
mulation The critical current density for spin- wave 
emission is reached when the difference in "chemical po- 
tentials" of the spin subbands equals the energy of a spin 
wave A/u = /ij — /ij = hu), so that in a spin-flip transi- 
tion of a conduction electron the energy is conserved [fl()| . 
In this way magnetic moment is transferred to the sys- 
tem of localized spins at the non-magnetic/ferromagnetic 
metal interface as the conduction electron spin polariza- 
tion changes. 

For spin wave emission a second ferromagnetic layer 
(e.g. the left layer in Fig. |l|) is not necessarily required 
as demonstrated in Fig. 2E of Ref. |l| . As long as only 
linear spin-waves are excited, the transfer of magnetic 
moment to the local moments leads just to a reduction 
of magnetization p6|] , similar to increasing the temper- 
ature. The reorientation of the axis of quantization of 
the ferromagnet, i.e. the switching, takes place only if 
non-linear spin-waves are created [|2l|| which occurs for 
much higher current densities. For currents below a crit- 
ical value the relaxation mechanism between the conduc- 
tion electrons and the local moments allows only for a 
damping of already existing spin- waves fic|| , similar as in 
Sec. V. 

Since switching seems to occur for current densities 
well below the critical value for spin- wave excitations, 
this mechanism is insufficient to lie at the origin of the 



switching. In particular for Co/Cu/Co spin valves, the 
data allows us to separate the hysteretic switching from 
the spin- wave excitations (refer to Fig. 2 in Ref. as 
mentioned earlier a spin-wave instability is observed al- 
ready for an initially unpolarized current. This figure also 
shows clearly that the domains are reversed well before 
the onset of spin-wave instabilities and at current densi- 
ties where transport yields more or less ohmic behaviour. 
Thus, experiments seem to provide evidence that the ex- 
citation of spin waves by the spin current is preceded by 
the switching process and, therefore, confirm our 
model for the switching based on the NEXI. 



VII. CONCLUSION 

In conclusion, we showed that the current driven cou- 
pling provides the dominant energy term that promotes 
the switching of the magnetization of the layers in spin- 
valves We pointed out the differences to the in- 
terpretation that the switching comes about by creating 
multiple spin-wave excitations due to inelastic spin-flip 
scattering which have to overcome only other forms of 
relaxation in the system. 

An interesting test of our model is to compare an asym- 
metric structure, such as Co/Cu/Py, with a symmetric 
one, such as Co/Cu/Py/Cu/Co. One could also sepa- 
rate the effect of spin- wave excitations, i.e. A/i = hio, for 
large currents from the switching, if one would apply an 
external field to the spin valve that is much stronger 
than /i^ exl given in Eq. (p^). On the other hand, inter- 
esting effects are expected to be observed in the region 
where h 1 ^* 1 and are comparable; in this case much 
of the magnetic response of the spin- valve is due to non- 
uniformity of the magnetization of the layers which we 
neglected here. 
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